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The existence of a geometric phase in magnetic traps can 
be used to Bose condense a magnetically trapped atomic gas 
into a vortex state. We propose an experimental setup where 
a magnetic trap together with a blue detuned laser beam 
form a multiply connected trap geometry. The local varia- 
tion of the magnetic quantization axis induces a geometric 
or Berry's phase that allows the atoms to acquire an effec- 
tive gauge charge interacting with the analog of a magnetic 
solenoid. It is shown that the ground state of such a system 
may be given by a vortex state. We also discuss the influ- 
ence of atomic interactions on the proposed vortex production 
scheme in the context of present Bose-Einstein condensation 
experiments with dilute gases. 
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The first experimental realization of Bose-Einstein con- 
densation in dilute alkali gases [|]-|| has dramatically re- 
vived the interest in weakly interacting Bose gases. Com- 
pared to previous work, the spatial inhomogeneity and 
finite size of the new experimental systems as well as 
the possibility to study complicated nonequilibrium dy- 
namics have added many interesting conceptual elements 
and have even shed new light on some basic paradigms 
of Bose-Einstein condensation. By now, also these novel 
Bose condensed systems seemed to be fairly well under- 
stood. However, as we will show in this paper, further 
surprises may still be abound. 

From the very beginning, the main focus of the investi- 
gation of Bose condensed alkali gases concentrated on the 
condensate itself, the macroscopically occupied ground 
state of the system. According to a credo of many text- 
books on quantum mechanics (cf. Q) the ground state of 
a particle in a spherically or axially symmetric potential 
has vanishing angular momentum. This statement is usu- 
ally derived from the famous theorem of nodelessness of 
the ground state || that in turn is based on the assump- 
tion that the eigenstates of the Hamiltonian can always 
be chosen to be real. However, this assumption does not 
necessarily hold for a charged particle in the presence of a 
magnetic field |f| . In our paper we show, that the action 
of an effective gauge field on neutral atoms in magnetic 
traps, caused by a geometric phase effect, can lead to a 
situation where the ground state of the system has a non- 
vanishing angular momentum and the gas can directly 
Bose condense into a vortex state. In so far, our vortex 
production scheme differs drastically from previous pro- 
posals that suggested the creation of a macroscopically 



occupied excited state of the system by opto-mechanical 
stirring [Q, adiabatic population transfer M or the acci- 
dental generation of vortices in a quench ]9|,|lfj[]. 

Before discussing the specific details of our proposed 
trap configuration, we want to illustrate in terms of clas- 
sical physics the basic mechanism that is used to cre- 
ate a ground state with nonvanishing angular momen- 
tum. Consider an electrically charged pendulum in the 
absence of gravity whose pivot point is pierced by an 
infinite solenoid perpendicular to the pendulum's plane 
of oscillation. The Hamiltonian for the pendulum then 
reads 
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where M and R are the mass and length of the pendulum 
respectively. / = MR 2 is the momentum of inertia of 
the pendulum, while C z = MR 2 <f> + AC Z is the angular 
momentum, canonically conjugated to the phase <j). The 
shift AC Z of the angular momentum is given by 



AL Z = 



q^B 
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where q is the pendulum charge, c is the speed of light, 
and $s is the magnetic flux across the solenoid. 

Apparently, the ground state of the Hamiltonian (|l|) 
posesses a nonvanishing canonical angular momentum 
C z = AC Z . In classical mechanics this increase of the 
canonical angular momentum does not lead to any im- 
mediately observable physical effect, since the angular 
frequency 4> = of the pendulum is the same as in the 
absence of a magnetic field. Only when the magnetic field 
within the solenoid is switched off, the conserved canon- 
ical angular momentum is transformed into mechanical 
angular momentum and the pendulum is set into mo- 
tion. Microscopically this generation of mechanical an- 
gular momentum is a consequence of the displacement 
current that is induced by a time-dependent magnetic 
field due to the Maxwell equation V x E = — (1/c)$b. 
In quantum mechanics, however, the canonical angular 
momentum is directly observable without altering the 
magnetic field. Here, the shift AC Z can be made visi- 
ble by an interference experiment similar to the one used 
to display the Aharonov-Bohm effect |ll[] . 

With this brief illustration we have shown how the ex- 
istence of a magnetic field can lead to a ground state with 
nonvanishing angular momentum. The direct use of the 
outlined idea in present experiments with Bose gases is 
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however flawed by the electric neutrality of atoms. In- 
stead, we show that the local variation of the hypcrfinc 
spin orientation in the inhomogeneous field of a magnetic 
trap leads to a geometric phase effect that can be inter- 
preted as the presence of an effective magnetic field. 

Consider an atom in a hyperfine state F, which moves 
in an external magnetic field B(r). The Hamiltonian of 
the system is given by 
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2M 

F(r)= 5MB P-B(r), 



(3) 
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with the atomic mass M, the Bohr magneton fis, and 
the Lande factor g. 

If the atomic motion is slow enough, the evolution of 
the system can be described using the Born-Oppcnhcimcr 
approximation. The latter assumes that for every point 
r the internal state of the atoms is close to a local eigen- 
state |m(r)) 



V(r)\m(r)) — m|m(r)) 
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of the magnetic contribution V"(r) to the Hamiltonian. 
The total wave function can then be written 



(r,m|*) = *^(r)(m|m(r)> 
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where |to) corresponds to an eigenstate of the projection 
operator F z — e 2 ■ F of the hyperfine spin on the z-axis. 

Applying the full Hamiltonian Eq. (|j) to the wave func- 
tion (r, m|^) of Eq. (^|) shows that the evolution of the 
adiabatic wave function is governed by the Born- 
Oppcnhcimcr Hamiltonian fll2fl 
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with the magnetic potential 

Vfn(r) = gfi B m\B(r)\, 
the vector gauge potential 

A J j l (r) = ih(m\Wm), 
and the scalar gauge potential 
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In addition to the phase change caused by the kinetic 
energy and the potential energy Vfh (r) the adiabatic wave 
function ^^(r) acquires the Berry's phase [Of 
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However, according to Eq. there exists an ambigu- 
ity in the definition of the adiabatic wave function ^m(r) 
and of the gauge potentials. The eigenstate equation (||) 
defines the local eigenstate \m(r)) only up to an arbitrary 
gauge transformation |m(r)) — > exp (r)] |m(r)). To 
eliminate this ambiguity we specify the gauge by choos- 
ing the following definition for |m(r)): 

|m(r)> = K d[^B(r)]| ' an §( e - B ( r ))) \ m = ™>' ( 12 ) 

The operator 7Z(n, 0) describes a three-dimensional rota- 
tion around an axis n by the angle 0. The state |m = to) 
is an eigenstate of the hyperfine spin projection on the 
z-axis with eigenvalue to. 

Our next task is to evaluate the potentials of Eqs. (||), 
(^), and ( |Io| ) for some particular magnetic trap. We 
choose a two-dimensional version of the Ioffe-Pritchard 
trap 



B(r) = Ax e x - Xy e y + B e z 



(13) 



assuming that the atomic motion is constrained to the 
xy-plane by some auxiliary external field that is indepen- 
dent of the internal hyperfine spin state. We will discuss 
a particular version of such a confinement below. We also 
restrict ourselves to the internal hyperfine state F = 1, 
to = — 1, g < that corresponds to present experiments 
with trapped sodium gases. The discussion of the gen- 
eral case is postponed to the end of the paper |L4j . Using 
polar coordinates for the iry-components of r, we obtain 
the expressions 
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on a semiclassical trajectory along the closed path T. 



with p = xe x + ye y and e ip (p) = (—ye x + xe y )/p. In 
the special case of a vanishing bias field Bq — 0, the an- 
gular momentum shift is equal to A£_i = —1 and the 
gauge vector field A_i(p) is equivalent to the electro- 
magnetic vector field of a solenoid along the z-axis with 
q&B/c = — 2-7T. In so far, the situation of the trapped 
atoms resembles the previous example of the electrically 
charged pendulum interacting with a solenoid. We there- 
fore conjecture that the ground state of the trapped atom 
gas has a nonvanishing canonical angular momentum. 
However, before coming to a definitive conclusion about 
the ground state of the system, we have to add a few 
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elements to the model outlined in Eq. (|13|) in order to 
make it realistic. 

The system Hamiltonian (^) does not yet contain a lon- 
gitudinal confinement along the z-direction. We therefore 
suggest to add a red-detuned light sheet parallel to the 
xy-plane. It has to be created by a far detuned and lin- 
early polarized light beam in order to achieve an ac Stark 
shift that is independent of the direction of the hyperfinc 
spin F. To a good approximation such a dipole trap can 
be described by a harmonic potential with oscillation fre- 
quency oj z . 

Since it is our plan to create a vortex in the absence of 
a bias field Bq, we have to find a way to suppress atom 
losses due to Majorana spin flips in the center of the trap. 
Therefore another, far blue-detuned light beam should 
be added that essentially plugs the region with the most 
critical trap losses (cf. (|]). This linearly polarized light 
beam is assumed to have a potential energy V p i ug (p) — 
Vpi ug ,oexp(-2p 2 /'(i;p lug ) with beam waist w p i ug . 

Finally, we have to account for interactions between 
the gas atoms. For dilute gases this can be done by 
adding the mean-field pressure 



= (ip\NH G p- 



(22) 



F mf (r) = 5 A^|vI/_ 1 (r)| 2 
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to the Hamiltonian (^). The coupling constant g — 
AttTi a/M is given by the s-wave scattering length a, 
while N is the total number of atoms. Furthermore we 
assume the axial confinement to be much stronger than 
the atomic interaction energy (Tilo z ^> <7 ^ „ 1 1 2 ) . In this 
case the atomic wave function 



*_ 1 (r)=V-i(p)-X (0) (2) 
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factorizes into a radial part ifi—i and the ground state 
of the axial potential x^°H z )- The radial wave func- 
tion ?/>_i(p) is then determined by a two-dimensional 
Gross-Pitaevskii Hamiltonian with an effective coupling 
strength g 2 v 



g\jMuT z j v 2-kTi. Note, however, that the 
existence of such an effectively two-dimensional situation 
is not instrumental for the proposed vortex creation. It 
was solely assumed to facilitate a numerical treatment of 
the problem. 

As a consequence, the ground state of the system 
is given as a solution of the two-dimensional Gross- 
Pitaevskii equation 



with the nonlinear Hamilton operator 
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the chemical potential p, the total atom number N. Our 
goal is to prove numerically that the minimum of the 
system energy 



is reached for a state with nontrivial value of the angular 
momentum C z = — 1. Using the separability of the radial 
and angular motion, we divide the energy minimization 
procedure into two steps: first, we minimize the energy 
under the constraint of a particular value of the angular 
momentum C z , then we compare the obtained energies 
in order to find its global minimum. 

As it was expected, we found that for the case of a van- 
ishing bias field Bq = 0, the energy minimum is reached 
for the C = — 1 state. As an illustration, we plot the 
energy difference between the C = — 1 and C — state in 
Fig. ffl. One can see that for a broad range of condensate 
populations the energy of the vortex state is constantly 
lower than the energy of the axially symmetric state. As 
a further illustration, Fig. ||a shows the radial distribu- 
tion of the atomic density of the vortex state. We regard 
this as a convincing proof of a possible Bose-Einstein con- 
densation into a vortex state. It is no surprise, however, 
that the described effect has not yet been observed ex- 
perimentally, since for typically used large bias fields of 
Bo S> Xp, p being a typical cloud radus, the condensation 
into a vortex state is energetically not favored. 

The main inconvenience of the proposed setup lies in 
the potential difficulties that it poses for the detection 
of a vortex state. The density distributions for different 
angular momenta are almost indistinct, since the most 
obvious signature of vorticity - a hole in the center of 
the density distribution - is masked by the optical plug. 
We therefore suggest the following experimental sequence 
in order to make the vorticity directly observable. 

• Turn on slowly a bias field Bq that is strong enough 
to prevent trap losses by Majorana spin flips. This 
stage is equivalent to turning off the solenoid field 
in the case of the pendulum Jl5| . 

• Remove the optical plug adiabatically, by lower- 
ing the light intensity (cf. Fig. |^b). Then, remove 
the magnetic confinement instantaneously and al- 
low the cloud to expand freely within the xy-plane. 

A detectable hole in the middle of the density distribution 
can now be taken as a distinctive feature of a vortex state. 
Alternatively, the vorticity of the state can be determined 
by an interference experiment fl6|| . 

In conlusion, we have pointed out that magnetically 
trapped atoms are subject to an effective gauge potential 
that is related to the existence of a geometric phase. Un- 
der certain circumstaces, the trapped atoms behave as 
if they were charged particles interacting with a magetic 
field of a flux carrying solenoid. Wc show that this effec- 
tive magnetic field can yield a ground state of the system 
with nonvanishing angular momentum, and thereby open 
the possibility of a direct Bose-Einstein condensation into 
a vortex state. Our results can be easily generalized to 
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an arbitrary value of the hyperfme spin F and to differ- 
ent kinds of magnetic traps. Indeed, one can show that 
the angular momentum of the ground state is then given 
by 



C z = 



(23) 



where T is the topological index of the used magnetic 
field configuration. It measures in terms of 2tt the ro- 
tation angle of the magnetic field vector B(r) along a 
closed path around the center of the trap. The proposed 
trap geometry of this paper corresponds to a topological 
index of T = — 1. In contrast, a quadrupole trap com- 
bined with an additional red-detuned light sheet in the 
xy-plane would correspond to an index T = +1 . 

After this work has been materially completed, we 
learned about a related vortex production scheme em- 
ploying the Ahoronov-Casher effect |L7| |. It combines a 
trap geometry very similar to ours with an electrically 
charged wire. In principle, such a scheme would allow 
the controlled creation of arbitrary vortex states, even 
though the experimental requirements seem to be chal- 
lenging. We show, in contrast, that Bose-Condensation 
into a vortex state may already occur in a more basic 
trap setup without the charged wire. 
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FIG. 1. In this figure, we plot for sodium atoms the energy 
difference AE = Ec z =o — Ec z =-i between energy minimiz- 
ing states with constrained angular momenta C z = and 
C z = —1 as a function of the atomic interaction strength 
(/2D-Y. We use the two-dimensional Ioffe-Pritchard trap with 
optical plug and confining light sheet that is described in the 
body of this paper. We use the realistic set of trap parame- 
ters A = 223Gs/cm, w = 3.05 /im (= 10 x {2h 2 //i B AAf ) 1/3 ), 
Upiug.o = 2ty x H x 9.7MHz, and u z = 26 kHz. For this set 
of parameters, the interaction strength g2oNM/?i 2 = 1000 
corresponds to N — 1250 atoms and an axial Thomas-Fermi 
parameter of (?2DW2D /hu) z pa 0.5, where fi2u is the peak value 
of the two-dimensional density in the xy-pl&ne. 
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FIG. 2. Here, we show the normalized two-dimensional 
density distribution n(p) for a vortex state with C z = — 1 and 
1250 sodium atoms (g 2D NM/h 2 = 1000). The left plot cor- 
responds to the interacting ground state of the toroidal trap 
described in the caption of Fig. [j]. The dot-dashed line shows 
the effective potential energy including magnetic field, light 
plug, gauge potentials, and the cenrifugal barrier. After the 
Bose-Einstein condensation into the vortex state is completed, 
a bias field Bo = 2 Gs is slowly turned on and the optical plug 
is slowly removed. The resulting trap configuration consists of 
a two-dimensional Ioffe-Pritchard trap with radial frequency 
uj±_ = 391 Hz, which corresponds to 0.08 in the units of Fig. 
1 ( and the axial light sheet described in the caption of Fig. 
1 The right figure demonstrates that the resulting atomic 
density distribution vanishes for small p as a consequence of 
the nonvanishing centrifugal barrier of a vortex state. 
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